Dissociative Electron Attachment to Polyatomic Molecules - I : 

Water 

N. Bhargava Ram,^'|jV. S. Prabhudesai,^ and E. Krishnakumar^'l^ 

^Tata Institute of Fundamental Research, Mumbai 400005, India 
^Department of Physics, Weizmann Institute of Science, Rehovot 76100, Israel 

Abstract 

Using the velocity map imaging technique, we studied and characterized the process of Dissocia- 
tive Electron Attachment (DEA) in polyatomic molecules like Water, Hydrogen Sulphide, Ammonia, 
Methane, Formic Acid and Propyl Amine. We present the details of these studies in a series of 5 
articles. In the first article here, we discuss the DEA process in gas phase water (H2O and D2O) 
molecules. Electrons of 6.5 eV, 8.5 eV and 12 eV are captured by water molecules in neutral ground 
state to form H2O (D2^ ) resonant states which dissociate into an anion fragment and one or 
I— I more neutrals. The details of the kinetic energy and angular distributions of the fragment anions H 

Qh (D ) and O produced from the three negative ion resonant states in the entire 27r scattering range 
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are discussed. Unique angular distribution patterns are observed at the 8.5 eV and 11.8 eV resonances 
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'^ showing dissociation dynamics beyond the axial recoil approximation. 
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I. INTRODUCTION 

Dissociative electron attachment (DEA) studies on water are of great significance for several 
reasons. It is the most basic molecule of life and hence, very important that we know the 
details of its interaction with radiation and charged particles that produce ions and radicals and 
consequent chemical reactions involving them. Recent studies have highlighted the contribution 
of the DEA process in damage to DNA and its significance in radiation induced damage to 
biological tissues [1-3J. Being the most ubiquitous molecule, water is also important as a 
solvent in many practical applications. Also being one of the simplest polyatomic molecules 
it is an ideal system to study and improve our understanding of the complex process of DEA 
in polyatomic molecules. Added to this, the availability of ab intio theoretical calculations on 
electron attachment to water pi-fTT] makes detailed experimental studies on the dynamics of the 
DEA process a timely necessity. The most important parameters in unravelling the dynamics 
of the DEA process are the kinetic energy and angular distribution of the products. Using the 
velocity map imaging technique, we measured these parameters for the H (D ) and O ions in 
the complete 27r angular range with unprecedented sensitivity. In this chapter, we present these 
data and identify the symmetry of the resonances, the dissociation channels and the dynamics 
at the three resonances. 

II. EARLIER WORK 

Electron attachment to water occurs as resonances centered at 6.5, 8.5 and 11.8 eV (Figure 
[I]) and has been widely studied to determine the exact resonance energies and widths, partial 
and absolute cross sections of the various anion fragments [12141^ . The electronic ground 
state of water is 1 ai^2ai^l b2 3ai^lb^^ -^ ^A^ state (02^ symmetry) and the three resonances 
are known to be core excited Feshbach resonances where an extra electron is attached to the 
singly excited states (caused by the excitation of the Ib^, Sa^ and 1 b2 electron respectively 
to 4aJ of the neutral water molecule. The short lived H2O dissociates to give H , O or 
OH fragments through various two-body and three-body channels as listed in Table 11} Also 
given are the possible symmetries of the H2O based on Wigner-Witmer correlation rules. 
The cross section for H ion production is maximum at 6.5 eV and lowest at 11.8 eV; whereas 
O" increases from 6.5 eV to 11.8 eV. The cross section for OH" ions is lower than O" by an 
order of magnitude at the above resonances [12j. The formation of OH has been a subject 
of controversy as some of the earlier measurements attributed its presence due to ion-molecule 
reaction [15j. Recent measurements by Fedor et al. [14j confirm its presence while theoretical 
calculations [H [TO] rule out the possibility of its formation by direct DEA. 
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FIG. 1: Cross sections for the formation of (a) H (D ) and (b) O from DEA to H2O (filled 
black squares) and D2O (open red circles) from ref. [13] 



TABLE I: Electron attachment to H2O and D2O : Possible dissociation channels with 
threshold energy and parent molecular states based on Wigner - Witmer correlation rules. 



From H2O 



From D2O 



Possible symmetry states 



H + 0H(X2n); 4.35 eV D + OD{X^U)] 4.52 eV A^, A2, B^ or B. 

H~ + 0H*(A2S); 8.38 eV D" + 0D*(A2S); 8.38 eV A^ or B2 

H + H + O; 8.75 eV D + D + O; 9.02 eV A2, B^ or B2 



O +H^; 3.56 eV 



O +D2; 3.68 eV 



Ai, Bi or B. 



O + H + H; 8.04 eV O + D + D; 8.28 eV 



Ai, Bi or B. 



OH (^S) + H; 3.27 eV 00(^2) + D; 3.44 eV 



Ai or B. 



The first measurement on angular distribution from DEA to water was reported by Trajmar 
and Hall |16]. They measured the angular distribution of H ions at 6.5 eV and identified the 
H2O " resonance as ^B^ decaying into H"(^S) + OH(X^H). The angular distribution was seen 
to peak at 100° with an asymmetry about 90° and dropped off very fast at lower and higher 
angles. Further, they suggested using symmetry arguments that the transition matrix element 
for H2O formation would be zero for electron beam axis in the molecular plane and non zero 
for axis perpendicular to the molecular plane and hence, the H-OH bond dissociates when the 
electron beam is perpendicular to the molecular plane leading to distribution around 90°. They 



could also infer the vibrational and rotational excitation of the OH fragment from the kinetic 
energy distribution of H . 

Belie et al. [17j measured the energy and angular distribution of H and D ions from H2O 
and D2O respectively using conventional 127° electrostatic analyzers. Based on their measure- 
ments, the three resonances were identified as having ^B^^, ^A^^ and ^^2 symmetry respectively. 
The energy distribution of the H ions at the first and second resonance revealed the diatomic 
OH fragment to be in the electronic ground state but with vibrational and rotational excitation. 
The rotational excitation of OH is shown to be more intense at the second resonance than in 
the first as the energy distribution of H" ions at 6.5 eV has better resolved vibrational structure 
whereas it is smoother in the second resonance at 8.5 eV. The third resonance was found to 
dissociate to the first excited state of OH(^I]+) and ground state H" ion. The electronically 
excited OH was also seen to be vibrationally and rotationally excited. 

The H" angular distribution data at the 6.5 eV by Belie et al. [17j agreed with the results of 
Trajmar and Hall [16] reaffirming the ^B^^ symmetry. Comparison with the electron energy loss 
(EEL) spectra of H2O measured by Chutjian et al. [IB] shows that the 1 b^^Ssa^^ -^ ^B^ excited 
state at 7.0 eV is the parent state of the DEA resonance at 6.5 eV. An extra electron in the 
3sa]^ results in a ^B^^ state. The shift in the angular distribution peak to 100° was explained 
as due to distortion of the incident electron wave from the plane wave due to direct potential 
scattering [T7] . Later theoretical work by Teillet-Billy and Gauyacq ^19j showed that the effect 
due to permanent dipole moment of the water molecule would explain the angular distribution 
exceedingly well. However, Haxton et al. [8] have recently shown that this backward shift is 
caused by a small contribution from the (i-wave to the main contribution from the p-wave and 
the subsequent interference between them. 

The angular distribution observed by Belie et al. ^ITj at the second resonance was said to 
correlate with ^A^^ state and result in H and OH fragments in ground state. The electronic 
configuration for this resonance was stated to be Sa^f^Ssa^ and said to match with the ^A^^ 
state situated at 9.3 eV in the EEL data [18]. Continuing on same lines, the third resonance 
was identified as second excited state of H20"^(^B2) plus two Ssa^ Rydberg electrons. The H" 
was accompanied with OH in the first electronic excited state (^S"^). In this way, the DEA 
resonances were identified as Feshbach resonances formed by the binding of two 3 sa^^ electrons 
to the first three electronic states of H20~^ formed by ejecting 1 b^, 3 a^ and 1 b2. This appeared 
to be consistent with the known strongly antibonding nature of the 3 sa^^ orbital making the 
respective potential energy surfaces repulsive leading to dissociation into H and OH fragments 
and with the calculations by Jungen et al. [20] . 

Curtis and Walker ^21j measured the kinetic energies of the D and O ions produced by 



DEA in D2O over the resonance energy regions of 6.5, 8.5 and 11.5 eV using a cylindrical 
mirror energy analyzer and a quadrupole mass spectrometer. The kinetic energy of D at 6.5 
eV showed the dissociation channel to be D"(^S) + OD(X^n) with vibrational excitation of 
OD fragment whereas the O was produced with near thermal energies and accompanied by 
D2 in highly vibrational excited state (O + D2; threshold 3.68 eV). They argued that D at 
the second resonance appears via the same dissociation channel as the first resonance but from 
a predissociation mechanism. It is known from experiments and models of photodissociation 
that the potential energy surfaces of the ^A^ and ^A^ states are alike. However, the ^A^ state 
of water correlates with H + OH(A^n) while the ground state OH(X^n), the major product in 
photodissociation arises from predissociation by the A ^B^^ state. Assuming a parallel correlation 
for the ^A^ anion [D" +OD(A^I])], they argued that the ^A^ anion state predissociates through 
the ^B^ state to give the ground state OD. According to them, O at the second resonance 
also arises from a sequential two body breakup leading to the O" + D + D (threshold 8.28 
eV) pathway. At the third resonance, two groups of D ions were seen; one coming from the 
D (^S) + OD(X^n) channel with kinetic energies ranging from 3 to 7 eV and the other group 
arising from a three body breakup (D + D + O) with energy partitioned in the ratio 0.4 : 0.4 
: 0.2. A similar sequential two body breakup was suggested to explain the production of O 
ions with the observed kinetic energies. 

In the last few years, there have been ab initio calculations of the water anion state potential 
energy surfaces by Haxton et al. [5HTTJ assuming ^B^^, ^A^^ and ^B2 symmetries. They deter- 
mined the various two body dissociation channels from the topology of the potential energy 
surfaces and calculated the cross sections for each of these channels taking into account the 
initial vibrational and rotational states of the neutral water molecule target. In the case of 
the first ^Bi resonance, they also calculated the angular dependence of the electron attachment 
cross section and determined the H angular distribution which agreed perfectly with the ex- 
isting measurements |16[ I17j . They could also reproduce the 100° peak and the asymmetry in 
the angular distribution in their calculation which they explain as mixing of partial waves. For 
the second and third resonances, they also take into account the Renner- Teller coupling effects 
(^A^ and ^B^) and the conical intersection (^B2 and ^A^) of the anion surfaces as they strongly 
determine the nuclear dynamics and consequently the cross sections. In one of their very recent 
reports [llj, they showed that the three body breakup is a major component of the observed 
cross sections in the O" channel at the second (^A^^) and third resonance (^B2) by subtracting 
the two-body breakup cross sections from the total cross section at the resonances. 

In the angular distribution calculations, they emphasized the role of the entrance channel 
amplitude in determining the electron attachment at a particular orientation [8|. The entrance 



amplitude for dissociative attachment is analogous to the dipole matrix element which controls 
the amplitude for photodissociation. The entrance amplitude, depends on two angles instead of 
one as in photodissociation. This function corresponds to the matrix element between a discrete 
resonance state and a background scattering wave function for an electron incident on the initial 
target state. The entrance amplitude depends upon the initial orientation of the molecule with 
respect to the incident electron beam and this dependence leads to an angular dependence in 
the cross section, even after it is averaged over the random orientations of the molecule with 
respect to the incident electron direction. Another factor is the "axial recoil" nature of the 
dissociation. Under, the axial recoil approximation, the orientation of the molecule in terms of 
the coordinates and (j) does not change as the dissociation occurs i.e., the recoil axis which 
connects the atom and the diatom center of mass does not rotate during the dissociation. If this 
is the case, then the probability for producing dissociative fragments at a certain orientation is 
the same as the probability for attachment at that orientation (if the survival probability for 
the dissociative species is unity) . If the approximation does not apply, then the dependence of 
the entrance amplitude upon the initial orientation of the molecule will effectively be spread 
over a range of final orientations, and final state- specific angular dependences are much less 
likely. The attachment probability is directly related to the laboratory-frame distribution when 
the axial recoil condition is met, requiring in the present context that the recoil axis which 
connects the atom and the diatom center of mass does not rotate during the dissociation. 
Using this procedure, the entrance amplitude has been shown to refiect the underlying shape 
of the Ib]^, Sa^^ and 1 b2 orbitals of neutral water, from which one electron is excited into 
the 4a]^ orbital to form the ^B^^, ^A^^ and ^63 Feshbach resonance configurations in a very 
recent paper [22j. This paper also describes experiments very similar to ours and reports 
kinetic energy and angular distribution measurements of H and O ions especially on the 
second resonance. Supported with calculations of angular distribution based on axial recoil 
approximation and classical trajectory simulations, it was shown that the resonance at 8.5 eV 
undergoes structural changes due to bending mode vibrations before dissociation thus distorting 
the angular distributions. 

Apart from electron attachment based studies, there have been studies/ measurements of 
the optical absorption and electron energy loss spectra of water that strongly facilitate under- 
standing DEA dynamics. Mota et al. [23j measured high resolution VUV photoabsorption 
spectra of water between 6 and 11 eV. Their data identifies the two absorption bands centered 
at 7.447 and 9.672 eV due to A^ -^ B^ and A^ -^ A^ transition respectively having mixed va- 
lence/Rydberg character and these excited states are known to be the parent states of DEA 
resonances in water at 6.5 and 8.5 eV. Electron impact measurements by Chutjian et al. |18] 



have revealed dipole allowed ^'^B^ (b^ -^ 3sa^) states centered at 7.0 and 7.4 eV (parent states 
of the ^B^ resonance anion), ' A^ (Sa^ -^ 3saJ states at 9.3 and 9.7 eV, optically forbidden 
transitions like ' A2 (b^ -^ 3pb2) at 8.9 and 9.1 eV and many others. A recent theoretical 
analysis by Rubio et al. [24j on the singlet and triplet excited states of water arising from 
the 1 b]^ and 3^^ molecular orbital, reveals that ' B^^ and ' A2 states are the lowest occurring 
electronic excited states with mixed valence/Rydberg character. 

III. ANGULAR DISTRIBUTION OF FRAGMENTS UNDER C2v SYMMETRY 

Water has C2V symmetry and the various symmetry states or representations (labelled A^^, A2, 
B]^ and B2) are defined by the four operations - Identity (I), rotation by 180° about the symmetry 
axis (C2), refiection in the molecular plane {a^) and refiection in the plane perpendicular to the 
molecular plane containing the symmetry axis (cr^). The character table for C2V point group 
and the basis functions written in terms of spherical harmonics transforming irreducibly under 
this group are given in Table [Tl| 

TABLE IL C2V group character and basis functions 



I C2 (Jv (y'v Basis Functions 



Al 1 1 1 1 1;" or y^^ + Yf"^] m even 



A2 1 1 -1 -1 Yl^ - Yf""; m even 



Bi 1-1-11 y^"^ + l^-"^; modd 



B2 1 -1 1 -1 Y^^ - Y-"^- m odd 



The method to calculate the angular distribution of fragments arising from a two-body like 
breakup for a polyatomic molecule of a particular point group has been described by Azria 
et al. [25j. We proceed on those lines to determine the angular distribution for H ions 
produced from dissociation of water anion of various symmetries under C2V group for partial 
waves (characterized by / values) whose capture is allowed. The quantity we are interested in 
calculating is | < Negativeionstate\Partialwave{\)\InitialState > p. Here, the initial state 
is the electronic ground state of water molecule i.e. A^ We make use of the Yq function to 
represent this ground state. The detailed steps leading to the angular distribution expression is 
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illustrated for the case of A^ (neutral) to B^ (negative ion) transition for the capture of p- wave 



(Yi). From Table II, the basis function for a B^ state is given by Yi + Y^^. Assuming that 
the dissociating H-OH bond is oriented at Euler angles (0, 6*, 0) with respect to the electron 
beam direction; we transform the partial wave Y-l^ from electron beam frame (or lab frame) 
along the dissociation axis. Appendix A lists the expressions for all Yf^s upto / = 3 rotated 
by Euler angles (0, 9, 0) leading to transformation from lab frame to dissociation frame. The 
expression arising out of the transformation of Y^ is given by equation (A. 3) in Appendix A. 
Similarly, we transform the basis function Y^ + Y-^^ defined for the molecular frame along the 
dissociation axis by Euler angles (0, /5, 0) using the equations (A. 2) and (A. 4) in Appendix A. 
Here, /3 is the angle between the molecular symmetry axis (C2 axis) and the dissociating 0-H 
bond. Assuming the ground state geometry of water molecule, /3 is half of H-O-H bond angle 
i.e. 52.5°. Rewriting the terms appropriately and simplifying the expression, the transition 
amplitude /(^,(/>) is given by: 



/(^, (p) (X < Y^ + Y^ ^\Y-I^\Yq > all terms in dissociation frame 

oc ^\/2sin6'sin0 (1) 



The scattering intensity lp{0) is given by: 



27r 

2 



1,(0) OC / |/|^# 



'0 



oc 2sin2 6' (2) 

Similarly, the expression for angular distribution due to a d-wave {Y2) capture under B^ 

state {Y2 + Y^^) is given by 

W) oc ^(0.64sin^^ + 0.36sin2 2^) (3) 

When more than one partial wave is captured in the process, interference between the indi- 
vidual amplitudes can lead to forward-backward asymmetry in the angular distributions which 
is generally the case. We obtain a general expression taking into account the interference terms 
with free parameters and use this expression to fit the experimental data. In the present ex- 
ample of B^ symmetry, the expression for the scattering intensity lpj^.d{0) taking both p and d 
partial waves is given by 



Ip+diO) oc 2a' sin^ d + 1.56^(0.64 sin^ 9 + 0.36 sin^ 26) + 2a6(1.05) sin sin 20 cos 5 (4) 

where a, 6 and 6 are the fitting parameters. The origin of 6 arises from the difference in the 
phases of the two partial waves due to potential scattering. 
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FIG. 2: Angular distribution curves for various symmetries under C2V point group for the 

lowest allowed partial waves. 



The angular curves for other symmetry states - A^^, A2 and B2 are calculated similarly using 
appropriate basis functions listed in Table [TT] and partial wave expressions and following the 
steps in the calculation detailed for the case of B^ above. The plots for individual lowest order 
partial waves for the four different symmetry states of the resonance when the neutral molecule 
is in Ai state, like in the case of water, are given in Figure [2j All the curves are normalized to 
their maximum value. The expressions for scattering intensity for each symmetry and allowed 
partial wave are listed in Appendix B. 



IV. EXPERIMENTAL SETUP 

The measurements were done using the new Velocity Map Imaging (VMI) spectrometer 
specially built to study low energy electron - molecule interactions, as described previously 
[26]. In this experiment, a magnetically coUimated and pulsed electron beam interacts with an 
effusive molecular beam formed by a capillary array. The negative ions formed in the interaction 
region are extracted into a time-of-flight mass spectrometer using a suitable extraction pulse 
after finite delay. The extraction field and the time-of-fiight tube along with a focusing electrode 
at the entrance of the fiight tube provide the needed VMI condition. The ions are detected 
using a two-dimensional position sensitive detector made of three 50 mm diameter micro-channel 
plates in Z-stack configuration and a Wedge and Strip anode. The ions striking the detector 
are recorded individually for their time of arrival (t) and position (x,y) in LIST mode using a 
CAM AC based data acquisition system running on LAMPS [27( l28j . The central slice of the 
Newton sphere which contains the velocity distribution of the ions with respect to the electron 
beam direction is obtained by selecting a narrow time window during analysis. 

For imaging H ions we had to make small changes from the previous set up in which we 
imaged mostly O" ions since for similar kinetic energies, H" ions travel 4 times the distance 
as O . In order to take this into account, we had to limit the electron gun pulse width as well 
as the delay between the two pulses. Also, the operating voltages for ion extraction, focusing 
electrode and fiight tube were increased appropriately for maintaining the VMI conditions. The 
increase in the bias voltage on the focusing electrode was found to affect the electron beam. This 
was eliminated by having a fine wire mesh on the puller electrode, which fianks the interaction 
region on the time-of-fiight side. The addition of this wire mesh modified the VMI conditions. 
SIMION simulations along with actual optimization restored the imaging quality. We also 
found that the magnetic field used for electron beam coUimation is affecting the imaging of H 
ions. Though we reduced the magnetic field strength to minimum level necessary for reasonable 
electron current, we find that the H ion images have some distortions due to it, which have 
not been corrected for. However, the D ions could be imaged without distortions. The kinetic 
energy of the ions is retrieved by calibrating the ion image size in terms of energy using the 
radius of O from DEA to NO or O2 at the same experimental conditions whose kinetic energy 
is known exactly. 

We are unable to separate the O and OH masses in our spectrometer. However, as the 
cross sections of OH at the resonances are lower than O almost by an order of magnitude, 
it is safe to assume that the angular distribution and kinetic energy distribution belong to O . 
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V. MEASUREMENTS, RESULTS AND DISCUSSION 

As discussed above, resonant electron attachment to water molecule appears to lead to the 
formation of transient H2O negative ion of B^^, A^^ and B2 symmetry respectively in the 
energy range 5-13 eV. These resonances are seen as peaks in the cross sections of the fragment 
anions produced from the dissociation of the water anion. As shown in Figure [T| the peaks for 
H" and D" ions occur at 6.5 eV, 8.5 eV and 11.8 eV respectively, while those for O" occur at 
7 eV, 9 eV and 12 eV respectively. The velocity images of H and O ions from H2O and D 
and O ions from D2O in the 5-13 eV energy range covering the three resonances are shown in 
Figures [3| [4| [5] and [6] respectively. The electron beam direction is from top to bottom through 
the centre of every image. The radial distance from the centre of the image represents the 
magnitude of the velocity of the ions when they are formed. Thus, the images give differential 
cross sections in both energy and angle. 

To obtain the kinetic energy distribution from the velocity image, ion counts are integrated 
over the 27r angular range as a function of radius and normalized thereafter. The radial coordi- 
nate is calibrated in terms of kinetic energy (eV) using the data of O" from DEA to O2 taken 
across its 6.5 eV resonance peak at the same experimental conditions as the fragment anion 
being imaged. O produced from DEA to O2 serves as a good calibration source as all the 
excess energy is carried away by the fragments as kinetic energy with no internal excitation. 
The expected kinetic energy of O for electron energies 5-8 eV ranges from 0.65 to 2.15 eV. The 
measured radius of the O images at these electron energies is converted into kinetic energy 
and then used to calibrate the radial coordinate to obtain the kinetic energy distribution of 
H and O . As the size of the velocity map image on the detector is independent of the mass 
of the ion that is imaged and depends only on the initial kinetic energy of the ions, the scale 
calibrated for the O ion cloud can be used for the hydride ions for the same voltage conditions. 
We also use the data of O" from DEA to NO for calibration especially when kinetic energies 
are less than 1 eV. The angular distribution is obtained by integrating ion counts over a range 
of the radius and in the angular range 9-5 to 9-\-b degrees covering scattering angles from 0° to 
360° in steps of 10°. 

A. First resonance process peaking at 6.5 eV 

The H" (and D") ions produced from the first resonance process (as seen in Figure pj- (a), (b) 
and (c) and Figure [SJ- (a), (b) and (c)) are ejected perpendicular to the electron beam direction 
with a continuous distribution of intensity along the radial direction which appear to increase 
with radius. The kinetic energy distribution of the H ions at 6.5 eV is plotted in Figure m^ a). 
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FIG. 3: Velocity images of H ions from DEA to H2O at various electron energies 




(a) O" at 6 eV 



(b) O" at 7eV 
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FIG. 4: Velocity images of O ions from DEA to H2O at various electron energies 

Figure [7|(b) shows H and D spectra in the same plot for comparison. Figure |7|[c) and (d) 
show the variation of KE distribution of H and D ions respectively at electron energies 5.5 
eV, 6.5 eV and 7.5 eV sweeping across the first resonance. 

As seen in Figure [71(a), the kinetic energy of H" ions range from to about 2.5 eV with 
a peak at 1.8 eV. For electron energy of 6.5 eV, the maximum kinetic energy observed is 2.5 
eV. Thus, about 4 eV energy is used in the dissociation process and the dissociation channel 
consistent with observed kinetic energies is H" + OH(^n) with threshold energy 4.35 eV. This 
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FIG. 5: Velocity images of D ions from DEA to D2O at various electron energies 




(a) O" at 6 eV 
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FIG. 6: Velocity images of O ions from DEA to D2O at various electron energies 



is in agreement with the observations of Belie et al. [17j. The broad KE distribution is an indi- 
cation of vibrational and rotational excitation of the OH fragment. In this context, one would 
expect discrete peaks in the kinetic energy distribution corresponding to these vibrational and 
rotational excitations. Belie et al. [TTj have reported observing clear peaks corresponding to 
the vibrational excitation of the OH fragment. However, we are unable to see distinctly the 
rings corresponding to these in the present VMI data due to the poorer energy resolution of the 
electron beam, which is about 0.5 eV. Though the energy resolution is not sufficient to clearly 
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FIG. 7: (a) Kinetic energy distribution of H" ions at 6.5 eV. The inner curves are the 

individual Gaussian fits obtained corresponding to vibrational excitations from v = 

(rightmost) to v = 4 (leftmost) of the OH radical, (b) Kinetic energy distribution of H" ions 

and D ions at 6.5 eV for comparison, (c) and (d) show variation KE distribution of H and 

D" ions across the first resonance process at electron energies 5.5 eV, 6.5 eV and 7.5 eV 

see the vibrational levels, we attempted to retrieve the intensity distribution of vibrational ex- 
citation by fitting for contributions corresponding to excitation of various energetically allowed 
vibrational levels of the OH fragment. The contributions of allowed vibrational excitations are 
taken as Gaussians with the same width, but different heights. The width and the heights are 
used as free parameters and the best fit for them is obtained. The fit shows occupation of 
vibrational levels i/=0 to v—^. The relative intensities of the vibrational levels obtained this 



way are given in Table HI along with the results of Belie et al. [17j and the calculations of 



Haxton et al. [6j. The present results show relatively larger intensities for the higher vibra- 
tional levels as compared to the previous measurement and the theoretical calculations. One 
reason for the difference with the previous measurement could be in the inadequate calibration 
of the transmission function of the energy analyzer used by Belie et al. [17j who optimized their 
spectrometer for relatively energetic ions around 1.7 eV. They mention that the transmission 
function of their kinetic energy analyzer was constant within 20% in the range of 1 to 3 eV. 
No information is given for the transmission function for their analyzer below 1 eV which is 



14 



expected to go down as with any conventional electrostatic analyzer. One of the advantages 
of a VMI spectrometer is its constant transmission throughout the energy range starting from 
zero, except for energies large enough to leave the image out of bound of the detector. The 
kinetic energy distribution of D ions is also similar to that of H ions except that it is scaled 
down due to the higher mass of D (by a factor 0.95) and slightly higher threshold of 4.52 eV 
as seen in Figure ^b). The KE distribution of H and D ions across the resonance as shown 
in Figures U{c) and (d) show the maximum kinetic energy increase with increase in electron 
energy from 5.5 eV to 7.5 eV. Also the distribution is broadened indicating excitation of higher 
vibrational states of OH (OD). 

TABLE III: Comparison of relative intensities of various vibrational levels of OH 



Vibrational Level Theory'^ Earlier Expt.^ Present Expt. 

Oil 1 

1 1.08 0.83 0.82 

2 0.60 0.30 0.57 

3 0.04 0.07 0.32 

4 - - 0.25 



"Haxton et al. [6j;^Belic et al. ^1 

The angular distribution of H ions at 6.5 eV over the entire kinetic energy distribution is 
shown in Figure IsFa). The intensity distribution for H appears to peak at an angle of about 
100° with respect to the electron beam direction. This is consistent with the data of Trajmar et 
al. [16j and Belie et al. [TTj as shown in Figure p[ a). From the observed angular distribution, 
we conclude the symmetry of the resonance to be B^^ in agreement with earlier reports. We 
fit our data with the functional form given in equation 3.4 for B^^ symmetry taking the lowest 
allowed partial waves - p and d. Resultant curve is shown in Figure |8| (a) as a solid line. We find 
the ratio of the d wave to p wave amplitude to be 0.31 with a phase difference of 3.139 (close 
to tt) radians. A similar procedure has been used by Haxton et al. [8j earlier to show that the 
observed asymmetry in the angular distribution at this resonance is due to the contribution 
from the d wave with a different phase. The angular distribution of D ions from D2O at this 
resonance is exactly similar to H ions. 
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FIG. 8: (a) H" angular distribution (integrated over all kinetic energies) at 6.5 eV (circles - 

present data, stars - Trajmar and Hall [1^, triangles - Belie et al.[l7j, solid line - fit for B^ 

symmetry using p and d partial waves). Angular distribution plots of (b) H^ and (c) D" ions 

as a function of kinetic energy at 6.5 eV. 



Figures p[b) and (c) show angular distributions of H and D ions plotted over smaller 
intervals of the kinetic energy distribution. All the plots peak close to 90° and fall off in 
intensity at lower and higher angles. We see no drastic change in the angular distributions 
except for a few observations. The plot for ions in the to 1 eV interval (shown by black dots) 
shows finite intensity at 0° and 180° compared to other plots. This is attributed to the fact 
that the ions with these kinetic energies are incident over a very small radial range (about 20 
pixels) compared to ions with 2 eV that hit at about 40 pixels distance. Considering energy 
resolution of about 0.5 eV of the electron beam and an extended interaction volume, poorer 
image resolution over this small radial range can lead to comparable intensities at forward- 
backward angles. Further, the angular distribution seems to be more symmetric about 90° at 
higher kinetic energies (or OH fragment in ground vibrational state) than at lower energies (OH 
fragment in higher vibrational states). 

The images for O" across the first resonance at electron energies 6eV, 7eV and 8 eV from 
H2O are shown in Figure |4|Ja), (b) and (c) and from D2O in Figure |6|^a), (b) and (c) respectively. 
They show a small blob suggesting very little kinetic energy which we are unable to resolve. 
The dissociation channel leading to the formation of O is O + H2 (threshold: 3.56 eV). The 
maximum kinetic energy expected for O is only l/9th of the available excess energy. If the 
H2 molecule is formed in the vibrational ground state, O will have energy close to 0.4 eV. In 
such a situation, we should be able to see a resolved ring for O as seen at the resonance at 
12 eV (as seen in Figures [Jg) and (h)). From the size of the image we estimate the kinetic 
energy of O" to be about 0.05 eV. This low energy implies that there is very little excess 
energy available as kinetic energy, with most of it going into vibrational excitation of the H2. 
This is easily understandable since the formation of H-H bond and the dissociation of the two 
0-H bonds are occurring simultaneously and most likely this will leave the H2 molecule with 
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large internuclear separation and hence large vibrational energy. In the case of 0~ from D2O, 
the expected maximum kinetic energy is 0.7 eV. We measure the kinetic energy to be 0.1 eV 
indicating that the D2 formed at this resonance will be in high vibrational states. Our results 
indicate that though low in cross section [13], the DEA to water around the 7 eV resonance 
produces D2 and H2 in very high vibrational levels. 

B. Second resonance process peaking at 8.5 eV 

The second resonance process peaking at 8.5 eV is understood to be a Feshbach resonance 
where an electron attaches to the (3 aj "^(4 aj ^ excited state of the neutral water molecule and 
forms a H2O " anion of ^A^^ symmetry. From the velocity images taken at 8.5 and 9.5 eV (see 
Figures|3|^d),|3je),|5jd) and|5]^e)), H" and D" ions are seen to be ejected mostly perpendicular 
to the electron beam (like in first resonance) with finite intensity in other angles as well. The 
kinetic energy distribution of H" ions at 8.5 eV is given in Figure ^[a). Here, the kinetic energy 
ranges from to above 4 eV. With incident electron energy of 8.5 eV, the only channel that 
allows maximum kinetic energy of H ions to about 4 eV is the H + OH(X^n) channel with 
threshold energy of 4.35 eV. This is the same dissociation channel producing H ions as at the 
first resonance. Like in the case of first resonance, we tried to retrieve the vibrational state 
intensities by fitting contributions at the allowed vibration levels. We find contributions of ten 
vibrational levels from i^ = to z^ = 9 are needed to fit the observed energy distribution. The 



relative intensities of these are given in Table IV along with those obtained by Belie et al [IT]. 



Though the maximum probability appears to be around i^ = 2 in both, the relative intensity 
distribution among the various levels appears to be different. The relative intensity distribution 
seems to go down faster in the data of Belie et al [17J with the increase in vibrational levels. 
The difference in the distributions could be attributed to the focusing of the energy analyser 
at 3 eV by Belie et al [17j. This may have led to discrimination of the lower energy ions which 
correspond to higher vibrational levels. The H kinetic energy spectra at 8.5 eV resonance as 
calculated by Haxton et al [lOj, shows the distribution peaking close to 3.5 eV which corresponds 
to ly = 1 vibrational state in the OH radical. The kinetic energy distribution we measured for 
D ions is similar to H ions at this resonance except for scaling down as seen in Figure [9Fb) 
due to difference in masses and slightly higher threshold of the D + OD(^H) channel i.e. 4.52 
eV. Looking at the variation of KE distribution, the plots at 8.5 eV and 9.5 eV for H and D 
ions in Figure pic) and p{d) show a slight increase in the maximum kinetic energy but there is 
no appreciable change in the width of the distribution. 

The angular distribution of H ions at 8.5 eV obtained by integrating over the entire kinetic 



energy distribution is shown in Figure 10 ^a) along with the previous results from Belie et al 
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FIG. 9: (a) Kinetic energy distribution of H ions at 8.5 eV. The inner curves are the 

individual Gaussian fits obtained corresponding to vibrational excitations from u — 

(rightmost) to z/ = 9 (leftmost) of the OH radical, (b) Kinetic energy distribution of D ions 

from D2O compared with that of H from H2O at 8.5 eV. Kinetic energy distribution of (c) 

H from H2O and (d) D from D2O at 8.5 eV and 9.5 eV electron energies. 



[TTj. The angular distribution that we have obtained has an isotropic contribution along with 
a bell-shaped distribution peaking at 90°. The data for D are very similar to that of H . As 



shown in Figure 10 'a), the angular distribution measured by Belie et al [17j peaks at 45° and 



135° respectively with finite intensity at 90° and corresponds to H" ions produced with OH in 
vibrational ground state (ly—O). They claimed that the peak at 90° due to H ions coming from 
the first resonance process and hence to minimise the contribution from the first resonance, they 
report the angular distribution of H ions produced with maximum kinetic energy. Based on 
these results, they concluded that the resonance to be ^A^^ state. The angular distribution from 
our velocity imaging measurement also peaks at 90°, but contribution from the first resonance 
arising due to the finite energy resolution of the electron beam (about 0.5 eV) was ruled out. 
The H angular distribution at this energy is distinctly different from that at the first resonance 
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TABLE IV: Relative intensities of OH vibrational levels at 8.5 eV 
Vibrational Level Theory" Earlier Expt.^ Present Expt. 

0.98 0.40 0.26 

1 L20 0.93 0.42 

2 LOO LOO LOO 

3 0.72 0.67 0.96 

4 0.51 0.41 0.88 

5 0.36 0.28 0.74 

6 0.26 0.24 0.46 

7 0.18 0.21 0.39 

8 0.12 0.14 0.20 

9 0.08 0.09 0.11 

"Haxton et al. [lOj; ^Belic et al. [El 

and measurement at 9.5 eV (where the contribution from the first resonance is negligible) is 
also found to be similar to that at 8.5 eV. 

Initially, we obtained the angular distribution from the velocity image by integrating over the 



entire radial range of the image (as shown in Figure 10 ^a)) and proceeded to explain our result 
in the following way. To determine the symmetry of the resonance, we fit our data for states of 
different symmetry, using the calculations for angular distribution as described in Section 3.3. 



The results for the A^^ symmetry are given in Figure 10 ^b). The three partial waves with lowest 
/-values which could contribute to the A^ state are the s, p and d waves. The probability that 
higher partial waves will contribute is increasingly small. It is clear from the figure that the fit 
obtained using these three partial waves are quite inadequate to explain our results. Fits for the 



B]^ symmetry using the p-wave and a combination of p and d waves are given in Figure 10 ^c). 
Of these two, the lone p-wave fits well the peak seen at 90°. However, this does not agree with 
the finite contribution that we see in the forward and backward directions. It might appear that 
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FIG. 10: (a) Comparison of H data at 8.5 eV from our setup and Belie et al. [17] (b) A^ and 
(e) B^ symmetry fits to our angular distribution data (d) Comparison with simulated angular 
distribution curves for the B^^ symmetry, taking into account rotational effects up to different 
angles - 0° (dotted curve), 65° (dashed line) and 90° (dash dot curve), (e) The red solid line is 
the fit with an s-wave contribution from a state of A^^ symmetry and a p-wave contribution 

from a state of B^^ symmetry. 



contribution from an s-wave should be able to explain this discrepancy. However, based on the 
symmetry, an s-wave contribution to a B^^ state is ruled out. Other possible symmetries like A2 
and B2 also do not explain the observed angular distribution. To summarize, the B^^ symmetry 
appeared to be the closest choice for explaining the results except for the finite contribution in 
the forward and backward directions. 

Next, we tried to account for this finite intensity in the forward-backward angles in the 
following way. The angular distribution fits discussed above is based on the axial recoil ap- 
proximation in which it is assumed that the dissociation takes place in a time scale before 
the molecule could undergo rotation or structural changes. The formation of the B^^ resonance 
requires that the electron approach the molecule normal to its plane. Thus, even if there are 
changes in the molecular geometry (due to bending mode vibration in this case) before the 
dissociation happens, the angular distribution should correspond to the p-wave cross sections 



shown in Figure 10 'c). However, if the lifetime of the state is long enough that the molecule 



undergoes rotation before it dissociates, then the angular distribution would get smeared out. 
We attempted a simulation of such a process by taking into account the effect of rotation before 
the molecules dissociate. We found that rotation about the principal axis by an angle of 65° 
would give a fairly good fit as seen in Figure floFd) . It may be noted that the bending mode 
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vibration coupled with rotation could smear out the distribution further. The rotation angle of 
65° in this case corresponds to a lifetime of 1.2 x 10~^^ sec for this resonance. Jungen et al [20 



calculated the lifetime of this state to be 6 x lO"^'^ sec, based on the A^ symmetry. That the 
resonance dissociates on a slightly longer timescale points to the fact that it may initially be 
formed in a bound potential energy surface and the dissociation occurs through coupling with 
other surfaces. We wish to point out that there is a slight asymmetry between 0° and 180° in 



the angular distribution data, which we could not account for. In Figure 10 'e), we also give a 
fit using contribution from both A^^ and B^^ states. Here, we consider only the lowest possible 
partial waves - 5-wave for A^^ and p-wave for B^^. The relative contributions from the A^^ and 
B]^ states were found to be 1:1.1. We find that this fit is in reasonable agreement with the 
data. Based on this, the presence of the A^ state cannot be ruled out at this energy. However, 
the Wigner-Witmer rule and the near-absence of OH at this resonance seem to rule out the 
contribution from an A^^ state. 

One of the puzzles in the DEA to water has been the very low cross section for the formation 
of O and OH as compared to that for H despite their lower thresholds and the larger electron 
affinities for O and OH as compared to H. The low cross section for the formation of O may 
be explained since it involves considerable structural changes including the breaking of two 
0-H bonds and the formation of the H-H bond at the lower energy resonances. However, 
the formation of OH is a simple two-body breakup and there appeared no reason for the 
observed low cross section. Fluendy and Walker [29] argued that while the unexpectedly low 
cross sections for the formation of OD (or OH ) at the 6.5 eV resonance could be ascribed 
to its B]^ symmetry preventing the formation of the OH anion in the electronic ground state 
based on Wigner- Witmer rule and that at the 11.8 eV resonance due to the dissociation of the 
OH anion, there is no valid explanation for the low cross section at the 8.5 eV resonance, if it 
has an A^^ symmetry. Hence, they felt that the observed low cross section for the formation of 
OD at the 8.5 eV resonance should make the assignment of the A^^ symmetry to it doubtful. 
Their argument that the presence of a state with the A^^ symmetry should lead to larger cross 
section for the formation of OH may be used to remove the ambiguity we observe in explaining 
the observed angular distribution in Figure [Io|[d). 

While the resonance at 8.5 eV is attributed to 3a;^ ^ 3sa^ excitation giving rise to the A^ 
symmetry, there is some evidence supporting the B^^ symmetry for this resonance. From the 
measurements of Chutjian et al [l8], we identify the resonance at 8.5 eV to be arising from the 
"^A2 (b^ -^ 3pb2) state. Their EEL data [18j on D2O at an electron residual energy of 0.5 eV 
and a scattering angle of 90° show a broad structure at 8.9/9.1 eV, which is identified as a dipole 
forbidden transition (b^ -^ 3pb2) causing the ' A2 state. With the extra electron captured in 
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the 3pb2 state, the overall symmetry would be B^^ since A2XB2 -^ B^^. Following the work 
by Curtis and Walker [HJ, Fluendy and Walker [29] did classical trajectory calculations to 
relate to the formers' results. Interestingly, they could explain all the essential features of the 
experimental findings [21j using only two anionic potential surfaces of symmetry ^B^^ and ^B2. 
They suggested that both the first and the second resonance in water be ascribed to the decay 
of a single ^B^^ state but with different partial wave contributions at the higher energy. 

Thus, based on these arguments, we concluded that the observed angular distribution is due 
to the B]^ state arising from a dipole forbidden transition in the neutral with the dissociation 
taking place at larger time scale so that rotational effects contribute [30] 

However, we realised that the above conclusion is not true. On careful analysis of the H 
and D velocity images, we found that the angular distribution varies with kinetic energy (or 
the radius of the image) . The angular plots as a function of kinetic energy are shown in Figure 



11 for 8.5 eV and 9.5 eV. We find that H" (D ) with higher kinetic energies (close to 4 eV; 
OH in ground vibrational state) peaks at 45° and 135° as observed by Belie et al [17] whereas 
ones with lower kinetic energies (OH in higher vibrational states) peak at 90° with sufficient 
intensity at forward and backward angles. Recent report by Adaniya et al. [22] also mentioned 
that the H2O molecular negative ion at 8.5 eV once formed undergoes structural changes 
due to bending mode vibrations before the dissociation could take place, thereby distorting the 
angular distributions. Their classical trajectory calculations using the potential energy surfaces 
constructed for '^A' state showed that the H-O-H bond angle opens up quickly following electron 
attachment and reorients it towards and beyond linear geometry at 90°. Trajectories in which 
one of the H atoms takes more of the available kinetic energy are found to undergo less bending 
motion before dissociation. Thus, what we learn is, the dynamics of H ions taking away 
most of the available excess energy as kinetic energy can be described well by the axial recoil 
approximation and the H ions with lower kinetic energies ejected at 90° as seen in the velocity 
images are due to water anion becoming linear due to bending mode vibrations and ejecting 
H ions in this geometry. 

On comparing our results with that of Adaniya et al. [22j , we find that their data for the 
kinetic energy release of 4 eV and the classical trajectory results are significantly different from 
the axial recoil model. In contrast, our results for 4 eV KER at 8.5 eV and 9.5 eV show peaks 
at 45° and 135° similar to what is expected for A^^ resonance under axial recoil approximation. 



as shown by the fairly matching fits in Figure 12 [a) - (d), taking the first three dominant partial 



waves (s, p and d). For smaller KER, the angular distribution tends to peak at 90°, as shown 



in Figure 12 [e). Here again, our results are more consistent with possible deviation from axial 



recoil approximation as seen by the better agreement with the classical trajectory calculation 
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FIG. 11: (Angular plots of H and D ions (normalized at 90°) from the second resonance 

process at 8.5 eV and 9.5 eV as a function of kinetic energy. These plots show variation in 

angular distribution as a function of kinetic energy suggesting structural changes of the water 

anion due to bending mode vibrations prior to dissociation. 



than the experimental results in 

The O formation from the second resonance process peaks at 9 eV as shown in Figure 
[IJb) and the images shown in Figures |4|[d) - (e) and |6|(d) - (e) are taken at 9 eV and 10 eV 
respectively. Like in the case of the first resonance, we are unable to see a resolved image for 
O in both H2O and D2O at this resonance too. The threshold for three-body breakup leading 
to O formation is 8.04 eV and earlier measurements [13j show that the entire O formation 
at this resonance occurs above this energy. The O + D + D channel was previously observed 
for D2O [21]. Thus we may conclude that the O is formed in this case entirely by a three- 
body fragmentation. Assuming instantaneous and symmetric three-body fragmentation, the 
O" kinetic energy would be equal to E'qCOS^ ^/(8 + cos^ 0) for H2O and E'qCOS^ 0/{A -\- cos^ 0) 
in the case of D2O where Eo is the excess energy and 9 is half the bond angle. The maximum 
KE would be £"0/9 and Eo/'o i.e. O.lleV and 0.19 eV respectively at incident electron energy 9 
eV for H2O and D2O respectively. For a bond angle of 105°, the energies work out to be 0.042 
eV and 0.081 eV respectively. From the velocity images, we estimate a kinetic energy of 0.06 
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FIG. 12: (a)-(d) A^ symmetry fits for H" and D" ions with 4 eV KER from the second 

resonance process at 8.5 eV and 9.5 eV. (e)H" ions with KE 2 eV at 8.5 eV electron energy 

and comparison with classical trajectory curve given by Adaniya et al. [i22j 



eV for O from H2O and 0.1 eV for O from D2O. The bond angles corresponding to these 
observed energies works out to be 86° in the case of H2O and 99° in the case of D2O. The 
smaller bond angles show that as the negative ion resonance dissociates, the two H (D) atoms 
are getting closer and that the dissociation limit may be the H2 (D2) vibrational continuum. 

Measurements by Adaniya et al. |22] show O from the ^A^^ resonance exhibiting an intense 
and broad distribution in the forward scattering angles and a less intense and narrow distri- 
bution in the backward angles for electron energies across the resonance. Our measurements 
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show a small unresolved blob. In the experiment by Adaniya et al. [22], the velocity images 
are obtained using the Abel inversion algorithm, whereas our measurements using the slicing 
technique are more direct. The O angular distribution reported by Adaniya et al. [22j peaks 
at about 30° and has a broad minimum about 110°, their theoretical prediction assuming ax- 
ial recoil approximation shows the peak at close to 150° and the minimum at 70°. Thus, in 
their case, the measured distribution and theoretical prediction are mirror images about the 
90° direction. They explain this distribution of O ions arising from the molecule scissoring 
backwards and ejecting the oxygen through the mouth of the H-O-H bond, regardless of the 
direction, while the hydrogen recoil in the opposite direction. Interestingly, this kind of a dis- 
tribution of O ions with forward scattering is a seen in the images of O produced from the 
third resonance process. It may be that the observed distribution of O" seen at energies 9.5 
eV and 10 eV in their data may be due to contribution from the third resonance process. It 
can be seen from the excitation function curves in Figure [lib) that these energies may have a 
contribution from the rising part of the third resonance structure peaking at 12 eV. 

C. Third resonance process at 11.8 eV 

The velocity map image of H (D ) fragment ions across the third resonance centered 
at 11.8 eV shown in Figure |3|^f), (g), (h) and Figure [SJJf), (g), (h) present a very striking 
behaviour. These images show H (D ) ions ejected mostly in the backward scattering angles. 
We notice three clear rings in the momentum images of H and D from H2O and D2O 
respectively. Energetically, these correspond to three different dissociation pathways of the 
molecular negative ion: 



7:^20"* ^H- + OH{X^U) (Channel 1) 

H- + OH* (A^S) (Channel 2) 

H- + H + (Channel 3) 



The kinetic energy distribution of H" (D ) ions (Figure 13 ^a)) shows the three dissociation 
limits as three distinct structures - the first one peaking at 1 eV, second one at 2.5 eV and 
third one appearing as a broad structure between 3.5 and 6 eV. The outermost ring in the 
velocity image (Figure |3|[g) and|5][g)) is faint, but its presence can be seen in the kinetic energy 
distribution as the broad structure from 3.5 to about 6 eV in the case of H from H2O and from 
3 to 4.5 eV in the case of D from D2O. The H + OH(X^n) channel which has a threshold 
of 4.35 eV will have an excess energy of 7.4 eV at 11.8 eV electron energy. The maximum KE 
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FIG. 13: (a) Kinetic energy distribution of H" and D" compared at 11.8 eV. (b) and (c) 
Variation of KE distribution of H and D ions respectively across the third resonance at 
electron energies 10.8 eV, 11.8 eV and 12.8 eV. (d) Kinetic energy distribution of 0~ from 
H2O and D2O at 12 eV. (e) and (f) Variation of KE distribution of O /H2O and O /D2O 

respectively at 11 eV, 12 eV and 13 eV. 

of H" would be 17/18 of the excess energy which is 7.0 eV in this case. The maximum KE 
from our measurements is 6.2 eV. This difference may be attributed to vibrational excitation 
of the OH fragment. In the case of D from D2O, the kinetic energy distribution appears to be 
considerably lower, even after taking into account the difference due to higher mass. The D 
intensity drops just above 4 eV and appear to trail off to 5 eV as compared to the maximum 
possible energy of 6.7 eV. The lower kinetic energy distribution indicates that the OD(X^n) 
states are formed at this resonance with populations in very high levels, even higher than the 
corresponding case of OH(X^n). 

The middle ring in the VMI data gives the peak at 2.5 eV in the kinetic energy distribution 
and corresponds to H ions produced with OH in electronic excited state (^S) with threshold 
energy being 8.4 eV. This distribution seems to extend up to a kinetic energy of 3.3 eV. The 
maximum kinetic energy of H in this channel will correspond to the formation of the OH (^S) 
in the vibrational ground state. This energy is 17/ 18th of the excess energy and works out to 
be 3.2 eV, which is in fair agreement with what we observe. The maximum of the H intensity 
distribution for this channel observed at 2.5 eV approximately corresponds to i^ = 2 vibrational 
excitation of the OH (^S) radical. The kinetic energy distribution of D from D2O peaks at 2 
eV corresponding to the OD (^S) channel and appears to extend approximately to 3 eV, which 



26 



is close to the expected maximum energy for this channel. We note that the low kinetic energy 
release in the case of both H2O and D2O in this channel indicates formation of the hydroxyl 
radical in high vibrational levels. This shift to higher vibrational levels appears to be more 
pronounced in the case of D2O as in Channel 1 and may be explained as due to the larger time 
needed for the dissociation of D2O negative ion state. 

The innermost ring in the H image corresponds to the kinetic energy distribution as a broad 
peak between 1 eV and 1.6 eV with the intensity falling off rapidly below 1 eV. For the three 
body breakup process (H" +H+0; threshold 8.75 eV), the excess energy would be 3.05 eV which 
would be distributed as translational energy amongst the three atomic fragments. Assuming 
instantaneous three-body fragmentation of the resonant state under axial recoil approximation, 
H" will have a kinetic energy equal to 4£'o/[8 + cos^ 0] where Eo is the excess energy and is 
half of the bond angle (105°, in this case). This works out to be 1.43 eV. If we assume that 
the molecule undergoes bending motion during the fragmentation the value of 9 will change. 
Assuming symmetric breaking of the two 0-H bonds, depending on the value of 9 the kinetic 
energy will be between 1.33 eV and 1.5 eV. We note that about half of the intensity distribution 
in the three body channel has energies below 1.3 eV. This effect seems to be more pronounced 
in the case of D from D2O. The D kinetic energy for the three - body channel peaks at 
0.7 eV, considerably lower than the minimum kinetic energy of 1.2 eV one expects from the 
instantaneous break up. 

One might account for the lower than expected kinetic energies of H" and D ions from the 
three-body break up channel as due to the formation of the O atom in the excited ^D state, is 
given by: 

H2O-* ^ H- +H + 0(^L')(threshold energy=10.72 eV) (Channel 4) 

The excess energy in this case will be 1.08 eV at incident electron energy of 11.8 eV. The 
energy range in which the H ions will appear in such a channel would be between 0.48 eV 
and 0.54 eV and for D the range would be between 0.43 eV and 0.54 eV. These energy ranges 
are much lower than the observed distribution in both H2O and D2O and hence we may rule 
out this channel. The only way to explain the lower kinetic energy distribution from the three- 
body break up channel is to assume a sequential fragmentation, instead of an instantaneous 
three-body break up. There are two ways this sequential fragmentation might occur: 
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H2O-* ^ H + OR-* ^ H + O + R- (Channel 5) 

^0 + H:^* ^O + H + R- (Channel 6) 

Calculations by Haxton et al. [9j indicate a possible decay mode for the B2 resonance through 
a H2" + 0(^D) channel. However, this would have higher energy threshold leaving less excess 
energy as in the case of Channel 4. Moreover, formation of H-H bond followed by its break 
appears too contrived and hence we may rule out Channel 6. Though we do not have much 
information on excited states of the hydroxyl anion, in order to explain the results, we need to 
assume that the sequential break up occurs through Channel 5. Curtis and Walker [21 J have 
pointed out that the three-body break up may be happening through the sequential process as 
in Channel 5. They attribute the formation of O also through such a process in which the 
extra electron is bound to the O atom. As to be discussed later, the analysis of the O" kinetic 
energy spectrum from D2O also shows contribution from sequential fragmentation via OD 
state. 

The kinetic energy distributions of H and D ions across the third resonance process are 
shown in Figure p^b) and (c). It is seen that the peak structure between 2.5 and 3.5 eV due 
to the OH(^S) (correspondingly OD(^I]) in D2O) diminishes in intensity while the intensity 
of the peak close to 1 eV increases. This indicates that the three body fragmentation channel 
becomes more favoured with increase in electron energy. 

The kinetic energy spectra of O from H2O and D2O respectively at 12 eV are given in Figure 



13;^d). As expected the kinetic energy of O from D2O is about twice that of O from H2O with 
peaks at 0.4 eV and 0.2 eV respectively. At the 12 eV resonance, the formation of O" occurs 
only through a three body fragmentation process which has a threshold of 8.04 eV. Since the 
excess energy available for this channel is 3.96 eV, we expect the O" kinetic energy distribution 
to extend up to 0.79 eV from D2O and 0.44 eV from H2O. The data seems to be in reasonable 



agreement with this. Figures 13 'e) and (f) show the KE distribution of O" across the resonance 
i.e at 11, 12 and 13 eV which show increase in maximum KE shifting towards higher values with 
increase in electron energy. Assuming symmetric fragmentation of the resonance leads to the 
formation of O" , its kinetic energy would be equal to Eo cos^ 6*7(4 + cos^ 0) in the case of D2O 
and Eo cos^ ^/(8 + cos^ 0) for H2O where Eo is the excess energy and is half the bond angle. 
For a bond angle of 105°, the kinetic energy works out to be 0.34 eV for O (D2O) and 0.18 
eV for 0"(H20). The peaks we see in the kinetic energy distributions are slightly larger and 
correspond to bond angles of 96° and 99° for D2O and H2O respectively. This indicates that as 
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the two 0-H bonds are breaking the two H atoms tend to get closer. One interesting thing we 



note in the kinetic energy distribution of O from D2O (Figure 13 ^d)) is the discernible drop at 
0.35 eV. In terms of bond angle, this corresponds to 105°. We explain this as a signature of the 
sequential fragmentation process seen in the kinetic energy spectra of H and D . As the O 
energy increases beyond 0.35 eV the bond angle decreases indicating that the potential energy 
surface of the resonance has its steepest descent towards larger 0-D distances and shorter D-D 
distances. This would mean that 105° would be a cut off angle beyond which no instantaneous 
three body fragmentation could occur as the system cannot climb up the surface. Hence the 
O intensity seen at energies below 0.35 eV has to be explained through a sequential process 
similar to what is seen in the case of H and D kinetic energy spectra. In this channel, the 
intermediate state is OD and this may decay through the formation of either D or O . The 
kinetic energy of the fragments coming through the sequential process would be smaller since 
the first D would have taken away a good fraction of the available excess energy. Hence the 
sequential dissociation may explain the observed O distribution below 0.35 eV. We note that 
the signature of the sequential process to be present in the O kinetic energy distribution from 
H2O as well. 

The absolute cross sections from H2O and D2O show isotope effects in H /D and O 
channels [13j. It is more pronounced in the O channel and present at all three resonances. In 
the H /D channel it is seen only at the third resonance. The cross sections for D formation 
is about 75% of that for H formation in this energy range [13]. We have made an attempt 
to quantify the isotope dependence in terms of the two-body to three body fragmentation 
channels. This is shown in Table |V] where we give the ratio of the two-body to three body 
fragmentation cross sections in the hydride ion channel at different electron energies across the 
third resonance. We note from Figure p^b) and (c) that at lower energy the two body channel 
is higher and decreases with respect to the three body channel as the energy is increased. We 
also note that the three body channel is relatively stronger in the case of D2O as compared to 
that of H2O indicating that isotope effect is more pronounced in the two body fragmentation 
channel. 

In addition to the appearance of three clear rings, the most striking aspect of the momentum 
distribution of the ions across the 11.8 eV resonance is the strong forward - backward asymmetry 
in both H /D and O ions. The angular distribution as seen from the images appears to be 
very unique with the H being ejected predominantly in the backward hemisphere. In contrast 
the O is ejected entirely in the forward hemisphere. The previous measurements by Belie et al. 
in which they identified only one kinetic energy component for H at this resonance appear 



to peak at 90° as well as toward 0°. They did not report any angular distribution data on O 
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TABLE V: The ration of two-body to three-body fragmentations in H2O and D2O in the 

hydride ion channel at the third resonance 



Electron energy (eV) H" (H2O) D" (D2O) 

10.8 1.97 1.39 

11.8 1.00 0.84 

12.8 0.68 0.71 



What we observe is considerably different from that reported by them. The results of Belie et al. 
pTj indicated the resonance to be a B2 state. We calculated individual transition amplitudes for 
p and d waves for a B2 resonance using the method described in Section 3.3 and calculated the 
angular distribution for H using these partial waves in an effort to fit the experimental results. 
We calculate the transition amplitude functions for A^ (neutral state) to B2 (negative ion state) 
for the lowest allowed p and d partial waves. Now combining the two individual functions along 
with the interference term, the general expression for scattering intensity /p_|_d(^) from a B2 state 
due to p and d wave capture is given by 2a^(2 sin^ /3 cos^ + cos^ (3 sin^ ^) + |^^(| sin^ 2/3 sin^ + 
cos^ 2/3 sin^ 20) + 2a6\/3(cos /3 cos 2/3 sin sin 20 — sin /3 sin 2(3cos0(3 cos^ ^ — 1)) cos 6 where a, b 
and 6 are the fitting parameters. This equation is fit to the angular distribution data and 
the goodness of fit determines the efficacy of the axial recoil approximation in descibing the 
dissociation process. The angular distribution for the three dissociation channels (both H and 



D ions) at 10.8eV, 11.8 eV and 12.8 eV are shown in Figure 14 The red curve in each plot 
depicts the best fit for the angular distribution data at 11.8 eV. The best fits at 10.8 eV and 
12.8 eV are very similar. The fits with p and d partial waves reproduce the forward-backward 
asymmetry in all three fragmentation channels. However, it is seen that except for the OH (^11) 
channel, we are unable to get a reasonably good fit. It is likely that all three channels may 
have effects due to deviation from axial recoil approximation. However, it is difficult to visualize 
how this deviation from axial recoil approximation will give rise to the strong forward backward 
anisotropy that we observe. We believe that the sequential fragmentation seen in the three-body 
break up channel may have smoothed the forward-backward anisotropy as compared to the other 
two channels. The angular distribution for O is shown in Figure p^d) along with a possible 
fit for a B2 resonance taking into account contribution from p and d partial waves. Assuming 
axial recoil, the O" ejects along the C2 axis (i.e. (3—0) and not along the 0-H bonds. Thus, 
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the functional form of the fitting expression is 2a^ sin^ 9 + |6^ sin^ 29 + 2a6\/3 sin 9 sin 29 cos 6. 
In this case also it is difficult to say that the fit is consistent with the experimental results, 
though we are able to reproduce the peak at about 60°. We find the ratio of the amplitudes of 
the pio d waves to close to 5 (inner ring), 4 (middle ring) and 2.5 (outer ring OH) for the three 
H channels with a phase difference close to n radians. The ratio of p to (i wave amplitudes 
for the O channel is close to 1.3 with a phase difference of about 1 radian. 

Since we are accessing the same resonance which is decaying through three different channels, 
the angular distributions should correspond to the same state with identical partial waves. In 
order to fit the forward-backward asymmetry, we have a phase difference of almost tt radians 
in the two partial waves for the H channels, while it is small for the O channel. The phase 
shifts among the partial waves are expected to be produced by the distortion of the incoming 
plane wave due to scattering by the molecule as a whole. One may argue that appropriately 
phase-shifted partial waves take part in the capture process leading to the observed angular 
distributions. This raises an interesting possibility that one set of partial waves with appropriate 
phase relation get captured to give rise to the H" ions, while another set with a different set 
of phases get captured at the same resonance, but leading to the formation of O ions. It is 
not clear if the complex dynamics of the nuclear and electronic motion leading to DEA could 
retain such fine details of the phase difference between partial waves in determining the final 
dissociation channel. 

It appears that the most probable way we get different angular distributions are due to 
deviations from the axial recoil approximation in the three different channels. The angular 
distribution for the OH(X^n) channel peaking at 180° with respect to the electron beam indi- 
cates that the electron is approaching along the H-0 bond and that mostly this particular 0-H 
bond is being broken. This is partly consistent with the fact that the probability for electron 
capture into the (lb2 4a^) B2 resonance will be maximum when the electron approach the 
molecule along the 0-H bond from the hydrogen side. In this orientation, we expect to see H 
to be ejected at 75° as well since either of the 0-H bonds could break following the electron 
capture. That we do not see a peak at 75° is intriguing. The only possible way this could 
happen is a strong asymmetric stretching introduced in the 0-H bond along which the electron 
is approaching the molecule. This asymmetric stretch, caused during the attachment process, 
makes the resonance dissociate along that particular bond, yielding H . 

The behaviour of the angular distributions in the two inner rings may be explained in terms 
of the structural changes that the molecular negative ion is undergoing after electron capture. 
For the outermost ring, the H" has the largest possible kinetic energy and we expect the 
dissociation to be the fastest, despite some of the energy going into vibrational excitation of 
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FIG. 14: Angular distribution of ions across the third resonance at three electron energies . 

(a) D from D + D + O channel (b) H from H + H + O channel (c) D from 
D~ + 0D(2S) channel (d) H from H + 0H(2S) channel (e) D from the D + OD{X^U) 
channel (f) H from the H + OH(X^n) channel (g) O from O + D + D channel and (h) 
0~ from 0~ + H + H channel. The red curve in each plot is the fit obtained for the 11.8 eV 

data using B2 symmetry functions. 
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the OH fragment as seen from the kinetic energy spectrum. Thus this channel will mostly 
follow the axial recoil approximation. As the three channels are assumed to be produced from 
the same resonant state, this distribution should also represent the anisotropy in the overall 
electron capture. For the middle ring corresponding to the OH(A^I]) channel, the H energy 
is considerably reduced to 2.5 eV due to the electronic and vibrational excitation of the OH. 
This reduction in energy corresponds to a slower dissociation process. In this timescale the 
vibrational motions of the molecule could force some of the fragmentation to take place along 
the 0-H bond that is away from the electron beam direction. Simultaneously, the bending 
mode vibration brings the two H atoms closer, producing the observed secondary peak at 80° 
and smearing out the abrupt drop seen in the angular distribution below 140° in the outermost 
ring. In the case of the three-body fragmentation channel, the bending mode vibration appears 
to smoothen out the anisotropy further. The sequential fragmentation present in this channel 
will also smooth the angular distribution. 

Although for the formation of O both the OH bonds need to be broken, one can deduce 
the angular distribution of O using the axial recoil approximation if the role of vibrations 
and rotations are neglected. Such comparison of the angular distribution is shown in Figure 



14 [^g) and (h). We may qualitatively explain the difference between the fit and the observed 
distribution in the case of O also as due to deviation from the axial recoil approximation, 
including sequential fragmentation. 

VI. SUMMARY 

1. First imaging experiment on study of DEA to Water. 

2. Velocity images of H"(D ) and O" ions obtained for incident electron energies 5-13 eV. 

3. First resonance centered at 6.5 eV has ^B^ symmetry state and agrees with previous 
measurements. 

4. Second resonance centered at 8.5 eV proceeds via ^A^^ symmetry and undergoes bending 
mode vibrations to linear geometry causing crossover to B^^ symmetry. Revealed by H 
angular distribution changing as a function of its kinetic energy (or OH vibrational state) . 

5. OH"(^I]) -\- H channel having the lowest threshold of 3.27 eV and favoured by a A^ 
resonance not seen. 

6. Third resonance centered at 12 eV has ^63 symmetry and shows strong anisotropy in the 
scattering of H and O ions. 
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7. H from the third resonance process is produced via three dissociation channels including 
a three body fragmentation channel and is ejected in backward angles. 

8. O produced via a three body breakup is scattered in the forward direction. 

9. Instantaneous and sequential fragmentation processes inferred from the kinetic energy 
spectrum leading to the formation of H and O . 

10. Deviation from axial recoil leading to complex dissociation dynamics observed in second 
and third resonance processes. 
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Appendix A: Spherical harmonics tranformed by Euler Angles (a, (3, 7) 

Spherical harmonics tranformed from one frame to another frame by rotation through Euler 
angles {a,(5,j) can be related as per the expression below (Ref: W Thompson, Angular Mo- 
mentum - An Illustrated Guide to Rotational Symmetries for Physical Systems, WILEY- VCH 
Verlag GmbH & Co. KGaA (2004)). The term on LHS i.e is the spherical harmonic in the new 
frame and the Yi^rn^ on the right are in the old frame; 



11,^.(^,0)= 5I^'-(^'«)^U'(«'/^'7) 



'^■m,m'\P) 



^{l + m')!(/ - m')\{l + m)\{l - m)\ 

(-l)^cos2^+^'-"^-2^(f)sin2^+"^-"^'(f) 



X 



E 



{I -\- m' — x)\{l -\- m — x)\x\[x + m — m')\ 



max{0, m' — m) < X < min{l — m,l -\- m') 

Following are the expressions for spherical harmonics upto / = 3 transformed by Euler angles 

(«,A0) 



Y2,-2{0e 



= Yr 



0,0 



^n,„ + e-"sin^(|)n,i 



e-!!^y,_,+eosW,o-e--^yu 



e" sin^(^)ri,-i + ^ri,o + e- cos\^)Y,,^ 



(Al) 
(A2) 

(A3) 

(A4) 



e'"" cos''(|)y2,-2 - e'" cos'(^) sin(/3)y2,-i + ^ sm^{0)Y2,o 



2\/2 



■e-" sin^(^) sm,3r2,i + e""° sin^(f )5'2.2 



(A5) 
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1^2,-i(^e,0e) = e^'"sin/5cos2(^)y2,-2 + e^"cos2(^)(cos2(^)-3sin2(^))y2,-i 



^2,o(^e,0e) = e 



B 



^2"^ sin2/5y2-2 + e^ 



(A6) 



2^2 
2^2 



sin 2/5^2,-1 + ^ l2,o 



2^2 



a/3 
sin 2/5^2,1 + e-^2"^ sin^ /5y2,2 



1^2,i(^e,<?^e) = e^'"sin/5sin2(^)y2,-2 + e^"sin2(^)(3cos2(^) -sin2(^))y2,-i 
■ ^ Sin2/5y2,o + e-^"cos2(^)(cos2(^)-3sin2(^))y2,i 



2^/2 
-e-^2"sin/5cos2(-)y2,2 



(A7) 



1^2,2(^e,</>e) 



1^3,-3(^e,0e) 



e^2"sin^(f )y2,-2 + e^"sin/5sin2(f )y2,-i + ^si^i^ ^Y^,^ 



2^2 



+e-^" sin /5 cos2(f )F2,i + e"^'" cos^(|)1^2,2 



(A8) 



(A9) 



J3a ^^^6//^- 



e^^"cos^(|)y3,-3 
-e^2«y6eos^(^) sin(^)y3,-2 + e^^TlS cos^(^) sin^l^)^,,-! 
-2^5 cos^(^) sin^(^)y3,o + e^ ^15 cos^l^) sin^(^)y3,i 
-e-^^"y6cos(^)sin^(^)y3,2 + e-^^"sin^(^)y3,3 
Y.MOeAe) = e^^"y6cos^(^)sin(^)y3,-3 + e^'"(cos^(f) -5cos^(^)sin^(^))y3,- 



(A10) 



2 
2 



2 



+e^"yi0cos^(-)sin(-)(2sin2(^) -cos2(^))y3,_i + ^sin^/^cos/^Kj^o 



2 
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+e-^"Vl0cos(^)sin3(^)(sin2(^)-2cos2(^))y3,i 



+e-^'"sin4(f )(5cos2(^) - sin2(^))y3,2 - e-^"y6cos(^) sin^(^)y3,3 (All) 



y^M^eAe) = e 



^^-/l5cos^(^)sin2(f)y3,-3 



+e^2«yi0cos^(^)sin(^)(cos2(^)-2sin2(^))y3,-2 



2^ ^2 

— sm - 

2^ ^2 

2'' ^2 



+e'"cos2(^)(cos^(^)-8cos2(^)sin2(^) + 6sin^(^))y3,-i 



sm 



-v^sin/?(cos'/?-— )y3,o 
+e-*°sm^(^)(6cos^(^) - 8sm^(^)cos^(^) + sin^(^))y3.i 
+e-'^°^cos(|)sm^(|)(sm^(^) - 2cos^(|))r3.2 
+e-*^"\/T5cos'(-) sin''(^)y3,3 



(A12) 
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1^3,0 (^. 



g^3« V g.^3 pY^_^ + e^2a V g.^2 ^ ^^g ^y^ _^ 

4 ' 2v2 

+e'"\/3sin/5(cos'/5 - ?^^)y3,-i + cos/5(l - - sin^ /5)y3,o 



sin^ (3 cos y^Fa 2 



-e-^"V^sin^(cos^ /5 - ^)y3 1 + e'^^"^ --^ 

4 ' ' 2a/2 

-e-^^"^sin3/5y3,3 



(A13) 



^3,l(^e,0e) = e^^-yi5cOS^(f)sin^(^)y3,-3 



+e^2"yi0sin^(^)cos(^)(2cos2(^)-sin2(^))y3,-2 
+e^"sin2(^)(sin^(f)-8cos2(^)sin'(|) + 6cos^(f))y3,-i 



+A/3sin/5(cos^/5 

,/5, 



sin^/5 



)^3,0 



4 

+e-^"cos2(^)(6sin^(^)-8sin2(^)cos2(^)+cos^(^))y3,i 



2 



2 



2 



-i2a 



10sm(^)cos^(^)(cos^(^)-2sin^(^))y3,2 



+e 



-•3°yi5cos^^)sm^(^)y3.3 



(A14) 



1^3.2(<^e,</-.) = e''"V^cos(^)sm^(^)y3,-3 + e'=^°sin^(^)(5cos^(^) - sm\'^))Y,.. 



2 
2 



+e'"V^sin3(-)cos(-)(2cos'(-)-sm'(-))y3,-i + ^sin"/3cos/3y3,o 



+e-'"\/10sm(^)cos"(^)(cos'(^) - 2sm2(-))F3,i 



+e-^°cos^(^)(cos^(^)-5sm^(^))r3,2 
-e-'^°V6sm(^)cos^(^)y3.3 



(A15) 



^-3,3(6., 't'e) = e'^"sm'(^)Y,,-3 + e•^"^/6sin^(^) 003(^)^3,-2 + e«^/T5cos^(^) sm''(^)y3,-i 



+e-'^"^/6 sin(|) cos^(|)K„2 + e"'^" cos<'(^)y3,3 



(A16) 
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Appendix B: Angular distribution curves for C2v point group 

Expression for angular distribution for various symmetries under C2V group taking lowest 
allowed partial waves. (5 is the angle between the dissociating bond and the molecular symmetry 
axis i.e. C2 axis. 9 is the scattering angle with respect to the electron beam direction in lab 
frame. 

Ai to Ai transition 



ls{0) = 1 (Bl) 

Ip{e) = sin2/5sin2^ + 2cos2/5cos2^ (B2) 

I^{e) = — (sinVsin^6' + sin2 2/5sin2 26') + i(3cos2/5-l)(3cos2^-l) (B3) 

16 2 

ls+p{0) = a^ + 6^(sin^/5sin^6' + 2cos^/5cos^6') +4a6cos/5cos6'cos(5 (B4) 

ls+p+d{0) = a^ + h^{sin^ (5siYi^ e + 2cos^ (5cos^ e) 

+c^(:;^(sin^/5sin^6' + sin2 2/5sin2 26') + i(3cos2/5-l)(3cos2 6'-l)) 
16 2 

+4a6 cos (5 cos cos 5i 

3 1 

+26c(- sin (5 sin 2(5 sin sin 20 ^ — cos (3{Z cos^ /5 — 1) cos ^(3 cos^ ^ — 1)) cos 82 

+ac(3 cos^ /5 - 1) (3 cos^ 9-1) cos((5i + (^2) (B5) 



A^ to Ao transition 



Id(e) = cos2/3sin^^ + sin2/3sin2 2^ (B6) 



Ai to Bl transition 



Ip{e) = 2sin2^ (B7) 

I^[e) = 1.5(sin2/5sin^^ + cos2/5sin2 2^) (B8) 

lp+d{0) = 2a2sin26' + 1.562(sin2/3sin^6' + cos^/3sin226') + 2a6y3cos/3sin^sin26'cos(5(B9) 
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Ai to B2 transition 



Ip{e) = 2(25m^/3cos^6' + cos^/3sin'6') (BIO) 

Id{e) = -(-sin2 2/5sin^6' + cos'2/5sin2 26'+^ ^{Scos^O-lf) (Bll) 

/p+d(6') = 2a^(25m2/5cos2^ + cos2/5sin2^) 

+ ^b'^{\ sin^ 2/5 sin^ ^ + cos^ 2/5 sin^ 2^) 

+2a6y3(cos /5 cos 2/5 sin 6' sin 29 - sin /5 sin 2/5co56'(3 cos^ 6* - 1)) cos (5 (B12) 
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